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Phase transitions in the imaginary chemical potential region are studied by the Polyakov loop 
extended Nambu-Jona-Lasinio (PNJL) model that possesses the extended Z3 symmetry. The 
extended-^ invariant quantities such as the partition function, the chiral condensate and the mod- 
ified Polyakov loop have the Roberge- Weiss (RW) periodicity. There appear four types of phase 
transitions; deconfinement, chiral, Polyakov-loop RW and chiral RW transitions. The orders of 
the chiral and deconfinement transitions depend on the presence or absence of current quark mass, 
but those of the Polyakov-loop RW and chiral RW transitions do not. The scalar- type eight-quark 
interaction newly added in the model makes the chiral transition line shift to the vicinity of the 
deconfinement transition line. 

PACS numbers: 11.30.Rd, 12.40.-y 



I. INTRODUCTION 



mat ion 



The lattice QCD (LQCD) simulations have become 
feasible for thermal systems at zero quark chemical po- 
tential (/i) [lj. As for /1 2 > 0, however, lattice QCD has 
the well-known sign problem, and the results are still far 
from perfection; for example, see Ref. 0] and references 
therein. 

Several approaches have been proposed to solve the 
sign problem. One of them is the use of imaginary chem- 
ical potential, since the fermionic determinant appear- 
ing in the Euclidean partition function is real in the 
case; for example, see Refs. 0, 0, [H IE 0] and references 
therein. If the physical quantity such as chiral condensate 
is known in the imaginary /x region, one can extrapolate 
it to the real \i region, until there appears a disconti- 
nuity. Furthermore, in principle, one can evaluate with 
the Fourier transformation the canonical partition func- 
tion with fixed quark number from the grand canonical 
partition function with imaginary /j, [H, [tj H3| • 

Roberge and Weiss (RW) [ll| found that the partition 
function of SU(Y) gauge theory with imaginary chemical 
potential [i — iO / ' (3 for fermion number 



Z(0) 



DyjDipDA^ exp 



(1) 



is a periodic function of 9 with a period 2it/N, that is 
Z{6 + 2irk/N) — Z(0) for any integer k, by showing that 
Z(6 + 2itk/N) is reduced to Z{9) with the Z N transfer- 
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where U(x,t) are elements of SU(iV) with U(x,f3) = 
exp(—2iTrk/N)U(x,0). Here ip is the fermion field, F^ v 
is the strength tensor of the gauge field A u , and /3 is the 
inverse of temperature T. The RW periodicity means 
that Z(9) is invariant under the extended Zjy transfor- 
mation 
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U^, A v -> UA„U- 
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Quantities invariant under the extended Zjy transforma- 
tion, such as the effective potential Q{9) and the chiral 
condensate, keep the RW periodicity. Meanwhile, the 
Polyakov loop <P is not invariant under the transforma- 
tion ([3]), since it is transformed as <P — > <pe~ l27Tk / N . In 
general, non-invariant quantities such as <P do not have 
the periodicity. This problem can be solved by introduc- 
ing the modified Polyakov loop \P(6) = •Pexp(i9) invari- 
ant under ([3]), as shown later. 

Roberge and Weiss also showed with perturbation that 
in the high T region dQ(6)/dO and <£((9) are discontinu- 
ous at 9 = (2k + l)ir/N, and also found with the strong 
coupled lattice theory that the discontinuity disappears 
in the low T region. This is called the Polyakov-loop RW 
phase transition in this P a pe r and it is observed in the 
lattice simulations || 0, H, H, 0] • 

Figure Q] shows a predicted phase diagram in the 9-T 
plane. The solid lines represent the Polyakov-loop RW 
phase transitions, and the dot-dashed lines do the chiral 
phase transitions predicted by the lattice simulation, al- 
though the result of the simulation is not conclusive since 
the current quark mass taken is much heavier than the 
realistic one 5 ~ 10 MeV. 

The Polyakov-loop RW transition is not an ordinary 
first-order deconfinement transition. Although both 
transitions are defined by discontinuities of <P, the lat- 
ter is a jump of the absolute value |<P| from almost zero 
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Fig. 1: The RW prediction on the QCD phase diagram in the 
0-T plane. The solid lines represent the Polyakov-loop RW 
phase transitions, and the dot-dashed ones correspond to the 
chiral phase transitions. 



to a finite value, but the former is a discontinuity of <P in 
its phase, as shown later. The discontinuity means that 
the Polyakov-loop RW transition is surely first-order in 
the phase. 

As shown later, the chiral condensate a is also not 
smooth on the Polyakov-loop RW transition lines (the 
solid lines of Fig. [T]) unless a is zero. This is called the 
chiral RW transition in this paper. This is also not an 
ordinary second-order chiral transition. The former is a 
jump of da / d9 from a finite value to its minus sign, while 
the latter is a divergence of dajdQ. 

As an approach complementary to first-principle lat- 
tice simulations, one can consider several effective mod- 
els. One of them is the Nambu-Jona-Lasinio (NJL) 
model [Hj]- Although the NJL model is a useful 
method for understanding the chiral symmetry break- 
ing, this model does not possess a confinement mecha- 
nism. As a reliable model that can treat both the chi- 
ral and the deconfinement phase transition, we can con- 
sider the Polyakov loop extended NJL (PNJL) model 

E3, HS> HB, Iio, E3, HM> H^> US, HM, 12^, US, HS> 

in which the deconfinement phase transition is described 
by the Polyakov loop. 

The PNJL model has already been applied to the real 
chemical potential region and many interesting results, 
in particular on the relation between the chiral and the 
deconfinement phase transition, are reported 
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the confinement mechanism shifts the critical endpoint 
of the chiral phase transition [U [H, Hi| to higher T and 
lower \i [2l|, HH • In order to confirm that these analy- 
ses are reliable, we should test the validity of the PNJL 
model by comparing the model results with the lattice 
ones. This is possible in the imaginary chemical poten- 
tial region where the lattice simulation is feasible. If the 
PNJL model is successful in reproducing lattice results in 
the imaginary chemical potential region, this will imply 
not only that the PNJL model is reliable for both the real 



and imaginary chemical potential regions but also that 
the lattice results can be reasonably extrapolated to the 
real chemical potential region by using the PNJL model. 

In the previous paper [30], we first applied the PNJL 
model to the imaginary chemical potential region and in- 
vestigated the phase diagram in the chiral limit. Among 
many effective models, the PNJL model has both the 
chiral symmetry and the extended Z3 symmetry needed 
to reproduce the RW periodicity. We showed with the 
extended Z3 symmetry that the Polyakov-loop RW tran- 
sition is first-order in the phase of & and the chiral 
RW transition is second-order. We also showed that 
both the crossover deconfinement and the second-order 
chiral transition take place in the 9-T plane, as ex- 
pected. In this paper, we make more extensive analy- 
ses on the four phase transitions by newly adding quark 
mass term and scalar-type eight-quark interaction to the 
PNJL Lagrangian and seeing their effects on the transi- 
tions. Throughout these analyses, we will find that re- 
sults of the PNJL model are consistent with all the lattice 
results. 

This paper is organized as follows. In §2, we describe 
the PNJL model. In §3, the extended Z3 symmetry, 
which plays a crucial role in our analyses, are introduced. 
In §4, numerical results are presented. Section 5 is de- 
voted to summary. 



II. THE PNJL MODEL 

The model we consider is the following two-flavor 
PNJL Lagrangian: 



C =q{i^ v D v - m )q 

+ G s [(qq) 2 + {qi lb fq) 2 ]-U{P[A]^[A\ 



,n (4) 



where q denotes the two-flavor quark field, mo does the 
current quark mass, and D v — d v — iA v — i^Sq. The 
field A v is defined as A v = 8^gA° a ^- with the gauge field 
A^, the Gell-Mann matrix A a and the gauge coupling g. 
In the NJL sector, r stands for the isospin matrix, and 
G s denotes the coupling constant of the scalar-type four- 
quark interaction. In §4, we will add a scalar-type eight- 
quark [IE Hil interaction to the PNJL Lagrangian to 
discuss the effect on the phase diagram. Effects of eight- 
quark interactions are discussed also in the three flavor 
NJL model [3l|. The Polyakov potential U, defined in 
(fT2")) . is a function of the Polyakov loop <P and its complex 
conjugate 



-TiL, 



with 



L(x) = Vexp 



*- = ^TrLt, 



dr A4 (x, r) 



(5) 



(6) 



where V is the path ordering and A± = iA . In the 
chiral limit (m = 0), the Lagrangian density has the 
exact SU(2) L x SU(2) R x U(l) v x SU(3) C symmetry. 
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The temporal component of the gauge field is diagonal 
in the flavor space, because the color and the flavor space 
are completely separated out in the present case. In the 
Polyakov gauge ,_L can be written in a diagonal form in 
the color space [15j : 



£ _ e i/9(<feA 3 +08As) 



diag(e 



(7) 



where <j) a = </> 3 + fis/VS, 4>b = —<}>3 + 4>a/V3 and 
4> c = -(0 a + 4>b) = -2^ 8 /V3. The Polyakov loop <P 
is an exact order parameter of the spontaneous Z3 sym- 
metry breaking in the pure gauge theory. Although the 
Z3 symmetry is not an exact one in the system with dy- 
namical quarks, it still seems to be a good indicator of 
the deconfinement phase transition. Therefore, we use # 
to define the deconfinement phase transition. 

Under the mean field approximation (MFA), the La- 
grangian density becomes 

£mfa = q{i^D^ - (m + S B ))q - U M (<r) - U{$, T), 

(8) 



where parameters are summarized in Table I. The 
Polyakov potential yields a deconfinement phase transi- 
tion at T = To in the pure gauge theory. In the previous 
paper [3(|, hence, To was taken to be 270 MeV predicted 
by the pure gauge lattice QCD calculation. However, the 
PN JL model with this value of To yields somewhat larger 
value of the transition temperature at zero density than 
that predicted by the full LQCD simulation [34], [35|, H|| . 
Therefore, we rescale Tq in §4. 
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TABLE I: Summary of the parameter set in the Polyakov 
sector used in Ref. [lg]- All parameters are dimensionless. 



The variables of <P, <P* and a satisfy the stationary 
conditions, 



where 



M 



G s a 2 . 



(9) 



Using the usual techniques, one can obtain the ther- 
modynamic potential per volume 



Q 



T\nZ 



V 

= -2N } 
1 



d 3 p 



(2n) 



3E{p) + -Tr c ln(l + Le~ pB ) 

P 



-Tr c ln(l + it e -^ + (p))j +U m ((t)+U(^$*,T), 

(10) 



dfi/d<P = 0, df2/d<P* = 0, dfljdo = 0. (14) 

The solutions of the stationary conditions do not give a 
global minimum fl necessarily; there is a possibility that 
they yield a local minimum or even a maximum. We then 
search a global minimum directly by varying the variables 
and check the solutions to satisfy (HU). The physical 
thermodynamic potential fl(6) at each 6 is obtained by 
inserting the solutions at each 9 into pip . 



III. EXTENDED Z 3 SYMMETRY 



where E(p) = y/p 2 +M 2 , E ± (p) = E(p) ± fj, = E{p) ± 
i9 / [3 and M = too + S s . In this paper, the thermo- 
dynamic potential per volume, fl, is simply called the 
thermodynamic potential. After some algebra, the ther- 
modynamic potential fl becomes [l8j 



fl 



2N 



d 3 p 



(2tt) 3 



3E(p) 



\ In [1 + 3(«P + Si* e -^-(p)) e -^-(P) + e -3^"(p)] 

P 

1 

U M +U 



■ In [1 + 3(<P* + <£ e -^ + (p)) e -^ + (P) + e -^ + (P)] 

(11) 



We use U of Ref. [18j that is fitted to a lattice QCD 
simulation in the pure gauge theory at finite T [32l . |33| : 



T 4 2 



6 V 



#3) 



i(#*<2>) 2 



Ht 



(12) 
(13) 



The thermodynamic potential fl of Eq. (TTTj) is not 
invariant under the Z3 transformation, 

<5(0) -» #(6»)e' i27rfe/3 , <£(6>)* -> <£(6») V 27rfc/3 , (15) 



although W of (fl2|) is invariant. Instead of the Z3 sym- 
metry, however, fl is invariant under the extended Z3 
transformation, 



<5(0)* #(0)*e i2 i\ 



(16) 



It is convenient to introduce the modified Polyakov loop 
if - = e l8 $ and = e~ l0 <P* invariant under the trans- 
formation (jTHJ) . The extended Z3 transformation is then 
rewritten into 



,±i0 , „±i6„±i 



w(ey -» #(0)*, 



(17) 
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and fl is also into 



Q 



2N f 



d 3 p 
(2tt) 3 



3£(p) + - In [1 + 3^e- 0B(p) 

P 



3^* e -2/3E(p) e /3MB + e -3/3E(p) e /3 M B] 

\ In [1 + 3<T e-^ (p ) + 3 ^ e -2^s(p) e -/3 Ml 



+ e -3/3E(p) e -/9 i u E 



-C/m 



& 2 (T)T 4 



4 a«) 

where /^b = 3/i = iW/[3 is the baryonic chemical poten- 
tial and the factor e ± ' 3MB is invariant under the transfor- 
mation (|17[) . Obviously, J? is invariant under the trans- 
formation (TTT|) . 

As a feature of the extended Z3 transformation (jTTJ) , 
the external parameter \x — i9 / (3 varies with (|17[) . and the 
modified Polyakov-loop is invariant under (|17[) and then 
not an order parameter of the extended Z3 symmetry. 
The same is seen for the chiral transformation in the 
chiral perturbation theory [3?], HI]. The Lagrangian of 
the theory is 



L, 



chPT 



1 

16' 



F 2 Tr 



+ ^Tr{M q (4 G + C/ NG )}+--- (19) 



with 



V: 



XG 



(20) 



where F 2 is a positive parameter, and —v, £ a and M q 
represent, the expectation value of the quark bilinear, 
the Nambu-Goldstone boson fields and the real matrix 
for quark masses in three flavor space, respectively. The 
Lagrangian (|19|) is not invariant under the chiral trans- 
formation 



XG 



exp i ^ 9a U NG exp 



a=l 



a=l J 



(21) 



since the second term including M q in the Lagrangian 
breaks the symmetry explicitly, where 9\ and 9^ are ar- 
bitrary real parameters of the transformation for the left- 
and right-handed quark fields. However, once the quark 
mass matrix M q is assumed to be transformed as 



exp Qa^a M q exp i ^ 9\X a , (22) 



0=1 



the Lagrangian (|19p becomes invariant under the chiral 
transformation, and v is not an order parameter of the 
extended chiral symmetry anymore. 



Although we have no order parameter for the extended 
Z3 symmetry, this symmetry leads many useful conclu- 
sions 30] as shown below. Under the transformation 
9^9 + 27rfc/3, (H]) keeps the same form, if &(9) and 
V(6)* are replaced by W{9 + 27rfc/3) and V{0 + 2iTk/3)*, 
respectively. This means that the stationary conditions 
for \P(9) and &(9)* agree with those for <P r (6» + 27rfc/3) and 
W(9 + 27rfc/3)*, respectively, and then that 

9wk 9irk 

#(0+ ) = W{9) and &(9 + )* = W(0)*. (23) 

3 3 

The potential fl depends on 9 through &(9), \P(9)*, 
a{9) and e 3ie . We then denote fl{9) by Q{9) = 
fl{W{9),^{9)* ,e 3ie ), where a{9) is suppressed since it is 
irrelevant to discussion shown below. The RW periodic- 
ity of fl is then shown as 



=fl{9), 



2nk 



)V 



(24) 



by using (|2"3"|) in the second equality. 

Equation (TT5)) keeps the same form under the trans- 
formation 9 — > —9, if &(9) and are replaced by 
*P(—9)* and W(—9), respectively. This indicates that 



&{-&) = &{9)* and <P{-9)* = &(0). 



(25) 



Furthermore, fl is a real function, as shown in (fT5)) . Us- 
ing these properties, one can show that 



fi{9) ={fl{9))* = fl(<P(9)*,<F(9), 



-3id 



=[2(&(-9),&(-9)*,e- il °) = Q(-9). (26) 

Thus, fl is a periodic even function of 9 with a period 
27r/3. The chiral condensate a(9) is also a periodic even 
function of 9, a(9) = <j(9 + 2nk/3) = cr(—9), because it 
is given by <r(9) = dfi{9) / dm®. Furthermore, the quark 
number density p v = —df2 / ' d{iT9) is pure imaginary and 
a periodic odd function of 9. 

The modified Polyakov loop W has a periodicity of (1231) . 
The real (imaginary) part of & is even (odd) under the 
interchange 9 <-> —9, because 



(#(0)+!p(0)*)/2 = Re[#(-0)], 
(&(9) - &(9)*)/(2i) = -Im[«?(-6 



Re[#(0)] = 
lm[&(9)] = 



where use has been made of (|25j) . Thus, the real 
(imaginary) part of W is a periodic even (odd) func- 
tion of 9. Similarly, the absolute value |<^| (phase <fi) 
of the Polyakov loop is a periodic even (odd) func- 
tion of 9, because \W\ = ^(Re^]) 2 + (Im[<f]) 2 (<j> = 
arctan(Im[t?]/Re[!f])). 

Since fl(9), &(0) and a (9) are periodic functions of 9 
with a period 27r/3, here we think a period < 9 < 2ir/3. 
In the region, periodic even functions such as fl{9), a (9), 
Re[<^(0)] and |<?| are symmetric with respect to a line 
9 = 7r/3. This indicates that such an even function has a 
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cusp at 9 = 7r/3, if the gradient dfl/dO is not zero. Such 
a cusp comes out in the high T region, as shown in §4 
with numerical calculations. This means that the chiral 
RW phase transition at 9 = ir/3 is second order. 

Meanwhile, lm[#(0)], <j> and p v are periodic odd func- 
tions of 9. This leads to the fact that these are discontin- 
uous at 9 = 7r/3, if the odd functions are not zero there. 
Thus, the Polyakov-loop RW phase transition appears 
as first order in lm[!^(0)] and <f> and as second order in 
Rc[<^(0)] and |t?|. The RW transition appearing in p v at 
9 = 7r/3 is also first order. These are seen in the high 
T region, as shown in §4. The orders of the Polyakov- 
loop and chiral RW phase transitions and the RW tran- 
sition appearing in p v are not affected by the existence 
of the current quark mass and multi-quark interactions 
introduced later, since the periodicity and the odd/even 
property of the physical quantities are not changed. 

The dynamical variables W and ^* are also invariant 
under the continuous phase transformation, 



field on lattice, respectively. In this case, the RW peri- 
odicity is expected to be naturally satisfied. In Table II, 
we summarize these symmetry properties in the three ef- 
fective models together with QCD and LQCD. Among 
the effective models, only the PNJL model has the same 
properties as QCD does. 



Theory 


Chiral 


1z 


extended 


Z3 RW periodicity 


QCD 


P 


PP 


P 


P 


LQCD 


P 


PP 


P 


P 


NJL 


P 


ND 


B 


B 


PNJL 


P 


PP 


P 


P 


3-d 3-state Potts 


B 


PP 


P 


P 



TABLE II: Summary of symmetry properties in several theo- 
ries. In QCD, LQCD, NJL and PNJL, we consider the chiral 
limit, mo = 0; P means "preserved", PP does "preserved in 
the pure gauge sector", B does "broken explicitly" and ND 
does "not defined". 



for an arbitrary real parameter a. 



$* (27) 

However, the factor 



6 ±/3mb(_ e ±3ie ) and the potential ft including the fac- 
tor are not invariant. If fl were invariant under (|27[) . 
the continuous symmetry would lead to a simple relation 
qr (9 + a) = &(9), that is <P{9 + a) = e- ia <P(9), that guar- 
antees that # is a smooth function of 9. When T is small 
under the condition that p is imaginary and W and \P* are 
not zero, the thermodynamic potential (|18[) is reduced to 



fl 



2N f 



dh 



(2tt) 3 



3£(p) + I In [1 + W e-^P)] 

P 



— In [1 + 3<f*e^ B(p) 

P 



+U M - 



(28) 



This has no explicit Pb{— 3i9T) dependence. Therefore, 
at low temperature, fl is approximately invariant under 
(f27|) and <P can rotate smoothly as 9 varies. At high tem- 
perature, however, effects of the explicit pb dependence 
are not negligible and <P can not rotate smoothly. Thus, it 
is obvious that the Polyakov-loop RW phase transition at 
high T is originated in the factor e ±/3AIB in (|18p . At high 
temperature, the continuous symmetry under the trans- 
formation (|27p is broken into a discrete symmetry, i.e., 
the extended Z3 symmetry, through the factor e ± @ m . 

It is easily seen that the ordinary NJL model respects 
the chiral symmetry but it does not preserve the extended 
Z3 symmetry. On the contrary, the 3-dimensional 3-state 
Potts model [Hi, S3, El El respects the extended Z 3 
symmetry and then has the RW periodicity, but it does 
not possess the chiral symmetry, since the model is a 
paradigm of QCD in the large quark mass limit. In lat- 
tice QCD (LQCD), to avoid the quadratic divergence, 
the chemical potential should be introduced just like the 
fourth component of an imaginary constant vector poten- 
tial [13, i.e., e^Ui or e'^ul, where a and U 4 (= e laAl ) 
are the lattice spacing and the fourth component of gauge 



Although the original NJL model does not has the ex- 
tended Z3 symmetry, it may be a good approximation to 
the PNJL model. In the mean field approximation, the 
fermionic part L7^ JL of the NJL thermodynamic poten- 
tial is given by 



f? f - 

"NJL — 



- 2N f N c 
1 



d 3 p 



(27T)* 



-^-( P ) ] + ^ ln [! 
P 



-0E+(p) 



(29) 



It is easily seen that ^p NJL (^ = <P* = I) = L?^ JL . There- 
fore, the fermionic part of the NJL thermodynamic po- 
tential coincides with that of the PNJL model in the per- 
fectly deconfinement phase. In the confinement phase, 
meanwhile, J? PNJL (<£ = <P* = 0) = L7^ JL is not held in 
general. Actually, it is also easily shown that 



4njl(* = = 0) = 



2m 
3 



T) 



2tt, 
3 



T) 



N 



(30) 



Therefore, in the confined phase, the fermionic part of 
the PNJL thermodynamic potential does not coincide 
with L?jvj JL but with the Z3 symmetrized average of L7^ JL . 
However, in the zero temperature limit, the thermody- 
namic potential of the PNJL model is reduced to that of 
the NJL model as 



^pnjl(T = 0) 



-QN f 



d 3 p 



(2tt) 3 

+ U u ((T,p v {T^0,n,a))+U(T 
= {2 mh (T = 0). 



E(p)-9(-E-(p))E-(p) 

0, <p,<r 



(31) 
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Therefore, the difference between the NJL model and 
the PNJL model is significant only in the intermediate 
temperature region. In both the high and the low tem- 
perature region, the two model give similar results for 
the physics concerning the chiral symmetry. It should be 
also remarked that the vacuum term in the PNJL ther- 
modynamic potential (fTTj) is the same as that of the NJL 
model. Therefore, the parameters included in the NJL 
sector of the PNJL model are the same as the param- 
eters in the original NJL model, if they are determined 
phenomenologically at T = jj, = 0. This also ensures 
that the original NJL model is a good approximation to 
the PNJL model. Thus, the difference between the two 
models is significant in the intermediate temperature re- 
gion. Actuall y, t he PNJL model makes the critical end- 
point [27], l28l |29|| shift to higher T and lower \x than the 
NJL model does HH, [Hi. 
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Fig. 2: T dependence of the chiral condensate normalized 
by cr(T = 0, fi = 0) and the absolute value of the modified 
Polyakov loop ^(9) at 9 = for three cases of To- Increasing 
(decreasing) functions of 9 denote |^| (a) for all the cases. 



IV. NUMERICAL RESULTS 

Since the NJL model is nonrenormalizable, it is then 
needed to introduce a cutoff in the momentum integra- 
tion. Here we take the three-dimensional momentum cut- 
off 

Hence, the present model has three parameters Too, A, 
G s in the NJL sector. Following Ref. [Ill, we use A = 
0.6315 GeV and G s = 5.498 GeV~ 2 , but vary m as a free 
parameter. If we put too = 5.5McV, these parameters 
reproduce the pion decay constant f n = 93.3MeV and 
the pion mass M w — 138MeV. 

Figure [3 shows the T dependence of the chiral conden- 
sate a and the absolute value of the modified Polyakov 
loop, at 9 = Q for three cases of To. The quantity 
(a) indicates that a crossover transition takes place 
at T D = 240MeV (T c = 261MeV) for T = 270MeV and 
at T D = 176MeV (T c = 221MeV) for T = 190MeV. 
We took the original value To = 270MeV in the previ- 
ous work [30j , but in this paper we take the rescaled one 
To = 190MeV. In this case the values of Tq and Tb are 
closer to 170~180MeV, the results of the two-flavor full 
LQCD simulation [34], [HI, [3(1, than those in the case of 
T = 270MeV, although the difference T c - T D becomes 
larger as T decreases, as shown in Fig. [2j 

Figure [3] shows the T dependence of the chiral conden- 
sate and the modified Polyakov loop for three cases of 
too. In the chiral limit too = 0, the chiral condensate 
shows that a second-order chiral phase transition takes 
place around 230 MeV not only for 9 — but also for 
9 =/= 0. The transition becomes crossover when too is fi- 
nite. Meanwhile, it is found from the modified Polyakov 
loop that the crossover deconfinement transition that 
appears around 180 MeV becomes sharper as Too in- 
creases. This behavior is consistent with the result of 
the 3-dimensional 3-state Potts model [3^, 0, Hi], \^ in 



which a first-order deconfinement transition is observed 
for large quark mass. 
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Fig. 3: The chiral condensate a normalized by a(T = 0, [1 = 
0) and the absolute value of !f'(6 l ) as a function of T for three 
cases of mo; (a) corresponds to 9 = and (b) to 9 — ^. 
Increasing (decreasing) functions of T denote \<P\ (a) for all 
the cases. 

Figure 2] shows the thermodynamic potential Q as a 
function of 9 in two cases of T = 170 and 200 MeV. 
The potential Q is smooth everywhere in the low T case, 
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but not at d = (2k + 1)tt/3 in the high T case. This 
result is consistent with the RW prediction and lat- 
tice simulation [3] on the 9 and the T dependence of the 
QCD thermodynamic potential. Qualitative features of 
fi at 9 = (2k + l)7r/3 are the same as in the chiral limit 
case 13011. 
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Fig. 4: Thermodynamic potential Q (in GeV 4 ) as a function 
of 8. The dashed line corresponds to the case of T = 170 MeV 
and the solid one to that of T = 200 MeV. 

Figure [5] shows the real and imaginary parts of ^(9) 
as a function of 9. In the case of T = 200 MeV, the 
imaginary part of &(9) is discontinuous at 9 = (2k + 
l)7r/3, while the real part of \P(9) is continuous but not 
smooth there. The transition appearing at 9 — (2k + 
l)7r/3 is the Polyakov loop RW phase transition. In the 
case of T = 170 MeV, meanwhile, both the real and the 
imaginary part are smooth everywhere. 

Figure [5] shows the absolute value |<?| and the phase <f> 
of the modified Polyakov loop as a function of 9. In the 
case of T = 200 MeV, the phase <fi(9) is discontinuous 
at 9 = (2k + l)7r/3, while the absolute value \\P{0)\ is 
continuous but not smooth there. In the case of T = 
170 MeV, meanwhile, both the absolute value and the 
phase are smooth everywhere. All the results on the 9 
and the T dependence of & are consistent with the results 
of lattice simulations @, 0, HI, IE 0- As an interesting 
feature, the transition appears as first order in the phase 
<t> and the imaginary part of X P(9) and as second order in 
the absolute value and the real part of &(9). Qualitative 
features of the Polyakov-loop RW transition are the same 
as in the chiral limit [30(. 

Figure [7] shows the chiral condensate a as a function 
of 9. In the case of T — 200 MeV, a has a cusp at 
each of lines 9 = (2k + 1)tt/3. Thus, the chiral phase 
transition of second order comes out at 9 = (2k + l)n/3. 
Meanwhile, in the case of T — 170MeV, there is no cusp 
at 9 = (2k + l)7r/3, indicating no chiral phase transition 
there. As a feature not seen in the chiral limit, a and 
the cusp do not vanish even in much higher temperature 
than T = 200 MeV. 

Figure [5] shows the imaginary part of quark number 
density p v as a function 9; note that the real part is 
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Fig. 5: The modified Polyakov loop ^(9) as a function of 9; (a) 
for the real part and (b) for the imaginary part. Definitions 
of lines are the same as in Fig. 4. 



always zero. In the case of T = 200 MeV, it is discon- 
tinuous at 9 — (2k + l)7r/3, indicating that the phase 
transition is first order. In the case of T = 200 MeV, p v 
is smooth everywhere. The 9 and the T dependence of j 
are consistent with the results of lattice simulations 
Comparing Fig. [5] with Fig. [5jb) , we see that the imagi- 
nary part of the quark number density has 9 dependence 
similar to the imaginary part of the modified Polyakov 
loop. This is natural, because the former is related to the 
fourth component of the vector current, while the latter 
to the fourth component of the vector field. 

Thus, the PNJL results are consistent with the lat- 
tice ones [3, S 0] , except that the temperature difference 
Tq — Td is considerably large in the former but small in 
the latter. Since the lattice simulations made in the imag- 
inary [l region have small lattice sizes, the results are not 
conclusive necessarily. However, it is important to check 
whether the PNJL model can reproduce the present lat- 
tice data in the finite imaginary fi region. This is dis- 
cussed below from a qualitative point of view. 

Our calculations have no free parameter, since the pa- 
rameters of the Polyakov-potential sector have been fixed 
to reproduce the results of the lattice QCD simulations 
and the other parameters in the NJL sector have been ad- 
justed to reproduce the empirical values of / ff and M n , 
i.e., f n = 93.3MeV and M„ = 138MeV. Therefore, we 
need a new parameter to improve this situation. It is 
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Fig. 6: The modified Polyakov loop \P(8) as a function of 8; 
(a) for the absolute value and (b) for the phase. Definitions 
of lines are the same as in Fig. 4. 



0.03 




Fig. 7: Chiral condensate a (in GeV 3 ) as a function of 
Definitions of lines are the same as in Fig. 4. 



reported that the value of the critical temperature is sen- 
sitive to the strength of scalar-type eight-quark interac- 
tion in both the NJL model HU and the PNJL one [25j [. 
Furthermore, the NJL and PNJL models with the eight- 
quark interaction can reproduce lighter cr-meson mass 
that may be phenomenologically more favorable [25], [2t| . 
Therefore, we can consider the scalar-type eight-quark in- 
teraction as a reasonable extension of the ordinary PNJL 
model. 



> 

u 
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CL 




Fig. 8: Imaginary part of the quark number density p v (in 
GeV 3 ) as a function of 9. Definitions of lines are the same as 
in Fig. 4. 



We add the eight-quark interaction (2f| H| 
G s8 [{qqf + (qi l5 rq) 2 ] 2 



(33) 



to Lagrangian The scalar part E s of the quark self- 
energy and the mcsonic potential Um are modified into 

E„ = -2G s a - AG s8 a 3 , U M = G s a 2 + 3G s8 a 4 . (34) 

Figure [5] shows the T dependence of the chiral conden- 
sate and the modified Polyakov loop with and with- 
out the eight-quark interaction. We set the parameter 
G s8 to 270.9 GeV" 8 to reproduce m a = 600 MeV at 
T = 9 = [291 ] . where m a is the er-meson mass. We 
see for 9 = that the eight-quark interaction makes the 
chiral transition temperature Tq shift down to the decon- 
finement transition temperature Td ~ 180 MeV. Similar 

effect is seen also for 9 = 5- 
b 

Figure ]W\ represents the phase diagram in the 9-T 
plane in the case of too — 0- The phase diagram is sym- 
metric with respect to each of lines 9 = kit/ '3 for any 
integer k. The dashed curve between D and E represents 
the deconfinemcnt transition of crossover, and the dot- 
dashed curve between C and F does the second-order 
chiral phase transition. For 9 =/= fc7r/3, thus, the chi- 
ral phase transition occurs at higher T than the decon- 
finement transition does. The solid vertical line start- 
ing from point E represents the Polyakov-loop RW phase 
transition. The chiral RW transition of second order oc- 
curs on the line between E and F. Above F, the chiral 
condensate is zero. Point F turns out to be a bifurca- 
tion of the chiral phase transition line, and point E is 
the endpoint of both the Polyakov-loop RW and chiral 
RW transitions. The chiral phase transition line between 
C and F is shifted down to the vicinity of the decon- 
finement transition line between D and E by adding the 
eight-quark interaction. In the case without (with) the 
eight-quark interaction, the temperatures at the point C, 
D, E, F are T c = 221(188) MeV, T D = 172(172) MeV, 
T E = 190(190) MeV and T F = 280(236) MeV, respec- 
tively. 
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Fig. 9: The chiral condensate a normalized by a(T = 0, [i = 
0) and the absolute value of ^(9) as a function of T with and 
without the eight-quark interaction; (a) for 9 — and (b) for 
= ^. Increasing (decreasing) functions of 9 denote \}P\ (a) 
for all the cases. 



Figure [TT] represents the phase diagram in the 6-T 
plane in the case of mo = 5.5 MeV. Again, the phase dia- 
gram is symmetric with respect to each of lines 9 = kn/3. 
The dashed curve between D and E represents the de- 
confinement transition of crossover, and the dot-dashed 
curve between C and F does the chiral transition of 
crossover. Also for this finite mo case, the chiral tran- 
sition occurs at higher T than the deconfinement transi- 
tion. The solid vertical line starting from point E rep- 
resents both the Polyakov-loop RW phase transition and 
the chiral RW phase transition of second order. Point E 
is the endpoint of both the RW transitions. As a view- 
point different from the mo = case [30(, point F is 
neither a bifurcation of the chiral phase transition line 
nor an endpoint of the chiral RW phase transition, since 
the chiral condensate a is always finite on the vertical line 
9 = 7r/3, that is, also on the line above point F. Just as 
in the case of mo = 0, the chiral transition line between 
C and F is shifted down to the vicinity of the decon- 
finement transition line between D and E by adding the 
eight-quark interaction. In the case without (with) the 
eight-quark interaction, the temperatures at the point C, 
D, E, F are T c = 218(194) MeV, T D = 176(176) MeV, 
T E = 190(190) MeV and T F = 282(244) MeV, respec- 
tively. 



Fig. 10: The phase diagram in the 9-T plane in the case 
of mo = 0; the eight-quark interaction is switched off in (a) 
and on in (b). The solid vertical line starting from point 
E represents the Polyakov-loop RW phase transition. The 
dashed curve between D and E represents the deconfinement 
transition of crossover, and the dot-dashed curve between C 
and F does the second-order chiral phase transition. The 
second-order chiral RW transition occurs on the line between 
E and F. 



Finally, we make a brief trial of the PNJL extrapola- 
tion, that is, the lattice data obtained in the imaginary 
fi region is extrapolated with PNJL to the real fi region. 
The absolute value of the Polaykov loop, \<P\, is measured 
for the two-flavor case in Ref. @, although the lattice 
simulation of Wilson dynamical quarks has a small size 
of 8 3 x 4. As the first step, we try to reproduce the lat- 
tice data [|| with PNJL. The screening mass of quark is 
estimated to be a few hundreds MeV in the simulation, 
so we simply change only mo to 100 MeV in PNJL cal- 
culations done above with A = 0.6315 GeV, G s = 5.498 
GeV" 2 , T = 190 MeV and G s8 = 0. Figure QJ presents 
|#| as a function of T in the case of 9/(ir/3) = 0.96, 
where the solid (dashed) curve corresponds to the re- 
sult of PNJL (LQCD). From a qualitative point of view, 
both the calculations give the same property that the 
Polykov-loop phase transition is crossover also for this 
case. From a quantitative point of view, however, LQCD 
gives a sharper transition than PNJL. 

Although more careful and extensive reproduction is 
necessary in future, now we do the PNJL calculation 
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Fig. 11: The phase diagram in the 0-T plane in the case of 
mo = 5.5 MeV; the eight-quark interaction is switched off in 
(a) and on in (b). The solid vertical line starting from point E 
represents both the Polyakov-loop RW phase transition and 
the chiral RW phase transition of second order. The dashed 
curve between D and E represents the deconfinement transi- 
tion of crossover and the dot-dashed curve between C and F 
does the chiral transition of crossover. 
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Fig. 12: The absolute value of <S> as a function of T plane in 
the case of 9/(tt/3) — 0.96. The solid curve is the result of 
PNJL with mo = 100 MeV, while the dashed one is the result 
of lattice QCD [fj. 



in both the real and imaginary fi regions with the pa- 
rameter set determined above. Figure [TST a) presents the 
chiral phase transition curves in the /i 2 -T plane. The 
solid (dot-dashed) curves represent that the phase tran- 
sitions are the first-order (crossover), while the dashed 
curve does the chiral RW phase transition of the second- 
order. The curve from point F to point A is smooth, 
as expected. For comparison, in Fig. 113( b) we present 
the result of mo = 5.5 MeV, where the other parameters 
are fixed. Thus, the phase diagram is largely affected by 
the change of tuq. This result encourages us to do more 
accurate LQCD similations with larger lattice sizes and 
small current quark mass and also to make more exten- 
sive comparison of PNJL calculations with LQCD ones. 
In the forthcoming paper, we will make such a PNJL ex- 
trapolation for the case of four flavor, since most LQCD 
simulations [4| are done in the case. 
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Fig. 13: The phase diagram in the p?-T plane for (a) 
m = 100 MeV and (b) mo = 5.5 MeV. The solid (dot- 
dashed) curves denote the first-order (crossover) chiral phase 
transitions, while the dashed curve does the chiral RW phase 
transiton of the second order. 

In the lattice simulations 0, [H, Q , the extrapolaton is 
made in a simple way in which the function 



T 



2ri 



(35) 



is assumed for the chiral transition curve and the coef- 
ficients a n are determined so as to reproduce the lattice 
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data in the imaginary \x region. In this paper the coeffi- 
cients a n are adjusted to the PNJL result of Fig. [TBTb) 
in the imaginary /i region. In Fig. 1141 four dashed curves 
labeled (l)-(4) represent results of the simple extrapola- 
tion for m = 1,2,3,4, respectively. The dashed curve 
of m — 4 still deviates from the PNJL curve in the real 
fi region. Thus, the PNJL curve includes higher-order 
terms than m = 5 that the simple extrapolation can not 
follow accurately. Furthermore, the simple extrapolation 
can not predict the position of the critical endpoint B. 
Thus, the PNJL extrapolaiton is superior to the simple 
extrapolation. 
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Fig. 14: The phase diagram in the fj?-T plane based on the 
extrapolation of p5|) . Dashed curves labeled (l)-(4) corre- 
spond to four cases of m = 1, 2, 3, 4 in (|35|) . respectively. 



V. SUMMARY 

In summary, we have studied the Polyakov loop ex- 
tended Nambu-Jona-Lasinio (PNJL) model with imagi- 
nary chemical potential. The phase diagram in the 9-T 
plane is studied in detail. Since the PNJL model pos- 
sesses an extended Z3 symmetry, quantities invariant un- 
der the symmetry such as the thermodynamic potential, 
the chiral condensate and the modified Polyakov loop, 
automatically have the Roberge- Weiss periodicity that 
QCD does. There appear four types of phase transi- 
tions; Polyakov-loop RW transition, chiral RW transi- 
tion of second order, deconfinemcnt and chiral transi- 
tions. The orders of the two RW transitions appearing 
at 9 — (2k + l)7r/3 are determined by the RW periodic- 
ity and the even/odd property of the modified Polyakov 
loop and the chiral condensate. As an interesting feature, 
the Polyakov-loop RW transition comes out as first or- 
der in the imaginary part and the phase of the modified 
Polyakov loop, but as second order in the real part and 



the absolute value of the modified Polyakov loop. These 
results are more informative than the RW prediction [llj 
and the results of lattice QCD @, 0, i, @, 0] . 

We have investigated effects of current quark mass uiq 
on the phase diagram. The orders of the ordinary chiral 
and deconfinement transitions depend on the presence or 
absence of current quark mass, as expected. In contrast, 
the orders of the Polyakov-loop RW and the chiral RW 
transition are not affected by the presence or absence 
of the current quark mass. The bifurcation point of the 
chiral transition line, which appears in the phase diagram 
in the chiral limit, disappears in the case of finite mo. We 
have also studied effects of the scalar-type eight-quark 
interaction on the phase diagram in the 9-T plane. It is 
found that the eight-quark interaction makes the chiral 
transition temperature shift down to the vicinity of the 
deconfinement transition temperature in the whole range 
of 9. 

The success of the PNJL model comes from the fact 
that the PNJL model has the extended Z3 symmetry, 
more precisely that the thermodynamic potential (|18p is 
a function only of variables, , e ±/3/iB and a, invari- 

ant under the extended Z3 symmetry. A reliable effective 
theory of QCD proposed in future is expected to have the 
same property in its thermodynamic potential. This may 
be a good guiding principle to elaborate an effective the- 
ory of QCD. 

Throughout all the present analyses, we can confirm 
that the results of the present model are consistent with 
the lattice results. In the analyses of the chiral limit, 
the present model is even more informative than the lat- 
tice simulation. In this paper, we have compared our 
results with lattice ones only qualitatively, since lattice 
simulations in the imaginary chemical potential region 
have been done mainly for the four-flavor case. In the 
forthcoming paper, we will make quantitative compari- 
son of the present model with lattice QCD. Determining 
the parameters of our model directly from the lattice re- 
sults in the imaginary chemical potential region, we will 
extrapolate the lattice results to the real chemical poten- 
tial region by using the present model. 
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